The helicity modulus in gauge field theories 
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We consider the 4d compact U(l) theory with Wilson action and characterize its phase diagram using the 
notion of electromagnetic flux, instead of the more usual magnetic monopole. Taking inspiration from the flux 
picture, we consider the helicity modulus (h.m.) for this theory, and show that it is an order parameter for 
the confinement deconfinement phase transition. We extend the definition of the h.m. to an Abelian projected 
Yang-Mills theory, and discuss its behavior in SU(2). 



1. Introduction 

Consider a theory characterized by the pres- 
ence of a phase transition as one parameter (3 (the 
inverse temperature or some coupling of the the- 
ory) varies through some particular value (3 C . A 
usual way to describe the transition is provided by 
the so-called order parameter, defined as a func- 
tion of [3 which is zero in one phase and non-zero 
in the other. The meaning of this definition is 
two-fold: on one hand, such a function has cer- 
tainly a point of non-analyticity, which reflects 
the underlying non-analyticity of the free energy 
of the system (in the thermodynamic limit); on 
the other hand, this behavior can sometimes be 
used to relate the phase transition with the spon- 
taneous breaking of some symmetry (implicitly or 
explicitly defined) of the model, thus providing an 
appealing physical picture of the transition itself. 

In this paper we introduce an order parame- 
ter for the 4d compact U(l) Abelian theory, and 
study it numerically for the Wilson action 

s = -pJ2 cos 9p w 
p 

where /3 is the inverse bare coupling and Op is the 
plaquette angle. This theory is characterized by 
a strong coupling confined phase and a weak cou- 
pling Coulomb phase; in the confined phase, in 
analogy with ordinary superconductors, the C/(l) 
gauge symmetry is supposed to be spontaneously 
broken. Our order parameter will make use ex- 
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clusively of the notion of electromagnetic flux (in- 
stead than the more usual magnetic monopole) 
which we now introduce. 

2. The flux in Abelian theories 

Let L be the lattice size of a Ad hypercubic lat- 
tice with periodic boundary conditions (p.b.c). 
Our definition of the flux $ M „ through a given 
((a, v) orientation is: 

^ = 75 E £[(0pUU* (2) 

(j»,»0 P^u 

where [9p]- ni7r is the plaquette angle reduced to 
the interval [—it, tt]. A double sum is present: 

• the internal ^2 P is the sum over the pla- 
quettes in a single plane, and is 2nk (k G Z) 
valued because of the p.b.c. . 

• the external average -p- planes' over au 
parallel planes of the given orientation, is 
non-trivial because the flux through differ- 
ent planes can change due to the presence of 
magnetic monopoles (herein lies the connec- 
tion with the usual picture). The allowed 
values for <I> M „ are thus multiples of 2tt/L 2 . 

Monitoring the flux distribution v(<f>) in the two 
phases, one observes that v{<j>) is Gaussian (cen- 
tered at <p — 0) in the confined phase, while it 
is peaked around multiples of 2it in the Coulomb 
phase (defining so-called flux sectors); in the ther- 
modynamic limit, tunnelling between flux sectors 
becomes completely suppressed. This behavior 
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Figure 1. Flux free energy vs. the external flux <p 
in the Coulomb phase (upper and middle curve, 
resp. (3 = 1.5,1.1) and in the confined phase 
(lower curve, (3 — 0.8). The fit using Eq.© is 
superimposed (/3r = 1.2, 0.71). The volume is 4 4 . 

already provides a characterization of the phase 
transition. 

3. Response to an external flux 

One can now ask what is the response of the 
system to an external electromagnetic flux, in the 
two phases. The corresponding flux free energy is 
a straightforward extension to the Abelian case of 
't Hooft's non- Abelian twist free energy Like 
the latter [!]J, it has a characteristic behavior in 
each phase. We impose an extra flux <f> G R to 
the following stack of plaquettes (this is only one 
possible choice): 

stack = {Op | fi = 1, v = 2; x = 1, y = 1}. (3) 
The partition function of the system becomes 

Z{(j>) = f D9e m ^^ cos{ep+ ^ + ^^ cosep) (4) 



where 'stack' is the complement of the stack, i.e., 
consists of all the other unchanged plaquettes. 
Alternatively (and more conveniently from a nu- 
merical point of view), it is possible to perform 
a change of variables such that the extra flux is 
spread through all the plaquettes with the given 
orientation (/i, is), leading to 

, f>(£ o,, „) co S (e P + ^)+X; — cos 9 P ) 

Z(<j>)= D6e planos (5) 



where '(//, v)' indicates all the other orientations, 
through which no extra flux is imposed. 

Z(<j>) is 2tt periodic in </>. FigQ] shows the be- 
havior of the free energy F(<fi) = — logZ(0) of 



the system in the two phases. In the confined 
phase (where the flux can freely change) the sys- 
tem is insensitive to the presence of an external 
flux (F((f>) = const.). In the Coulomb phase the 
free energy is described perfectly by the ansatz 

F(0) = -log5>-^- 2 ^ 2 (6) 

k 

which says that around each flux-sector the de- 
pendence on the external flux is quadratic, as ex- 
pected classically, and that we must consider the 
contribution of all the sectors at the same time. 
Pn is the only parameter in this equation. It re- 
places (5 in the classical expression © , and there- 
fore plays the role of a renormalized coupling . 

4. The helicity modulus 

The last step in the construction of the order 
parameter is to note that the physical informa- 
tion contained in Fig^can be rephrased in a more 
concise way: if instead of the whole curve (as a 
function of <fi) we consider only the curvature of 
F{4>) at the origin (or at any other point), we get 
a function of [3 which is always zero in the con- 
fined phase (F(4>) = const.) and is different from 
zero in the Coulomb phase; that is, it is an order 
parameter. This construction was already known 
in the context of the 2d XY model 0] , where the 
name 'helicity modulus' was first introduced. In 
our context we define the helicity modulus 

which can be related to /3r via Eq.©. Comput- 
ing explicitly the double derivative (with F(<f>) 
defined as per Eq.©), one gets 

planes planes 

In Fig|2] we plot this observable for different 
volumes: its behavior is qualitatively consistent 
with an order parameter. Moreover, it is shown 
in P| that the drop to zero in the confined phase 
is exponential in (/3 C — /3), providing convincing 
evidence that the transition is 1 st order. 

5. Extension to non- Abelian theories 

We now turn our attention to Yang-Mills theo- 
ries, which display a finite temperature transition 
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Figure 2. The helicity modulus vs. (3. 

(say, at T c ) between a confining and a deconfincd 
phase; here we consider the gauge group SU(2). 

If we assume that (I), following [S], the long 
range properties of the theory can be well de- 
scribed by an effective Abelian theory, then we 
expect that the Abelian gauge ensemble obtained 
after a suitable gauge fixing and Abelian projec- 
tion must change from confining to deconfincd 
across T c . Therefore, the helicity modulus mea- 
sured in the projected ensemble should present a 
discontinuity at T c . Moreover, under the further 
assumption (2) that the effective action which de- 
scribes the projected ensemble is Wilson-like 

S e g = /3 c ff E cos 6» Pproj . (9) 
p 

where f3 e g is an effective coupling whose value 
could be determined by Inverse Monte Carlo, the 
helicity modulus is the same as in Eq.l(8j. which 
we rewrite for convenience as 

h(/3) = - P 2 CS H 2 {P) (10) 

H M = ^{ E cos ^> ( n ) 

(M, ") 
planes 

w) = ^(( E sin ^) 2 )- ( 12 ) 

Cm, v) 

planes 

To perform the Abelian projection, we use the 
Maximal Abelian Gauge 0, defined by 

<j 3 ) maximum (13) 

where the maximum must be found over the set of 
all gauge transformations. Not knowing the value 
of Peg, we measure separately H\{($) and 7?2(/3). 
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Figure 3. The helicity modulus Ea. (|I0|) (non- 
Abelian case) for a test value f3 e g = 1.0. 

In fact, we can try to infer (3 c g from the require- 
ment that Eq. (|I0|l represents an order parameter. 
The value (3 c g = 1.0 around the transition region 
works well, as shown FigOU 

The similarity with Fig is remarkable, but fi- 
nite size effects here are more suggestive of a sec- 
ond order phase transition, as they should be. 

This behavior can be considered as a strong 
support of both assumptions (1) and (2). Our 
measurements of the same observable Eq.(ll) af- 
ter fixing to another gauge do not show a sim- 
ilar behavior. This presumably singles out the 
Maximal Abelian Gauge as yielding an effective 
action particularly local and close to the Wilson 
action. Fixing to another gauge, the effective ac- 
tion contains sizeable additional terms, and the 
expression for the helicity modulus Eq.Q differs 
appreciably from Eq. l|lU|) . 
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